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Abstract. This paper is concerned with the incompressible limit of the compressible magnc- 
tohydrodynamic equations with periodic boundary conditions. It is rigorously shown that the 
weak solutions of the compressible magnetohydrodynamic equations converge to the strong so- 
lution of the viscous or inviscid incompressible magnetohydrodynamic equations as long as the 
latter exists both for the well-prepared initial data and general initial data. Furthermore, the 
convergence rates are also obtained in the case of the well-prepared initial data. 



1. Introduction 

Magnetohydrodynamics (MHD) studies the dynamics of compressible quasineutrally ionized 
fluids under the influence of electromagnetic fields. The applications of magnetohydrodynamics 
cover a very wide range of physical objects, from liquid metals to cosmic plasmas. The compressible 
viscous MHD equations in the isentropic case take the form (see, e.g., [HllTSll^T] ) 

dtp + div(pu) = 0, (1.1) 

dt(pu) + div(pu <8> u) + VP = (curlH) x H + /xAu+ (£ + A)V(divu), (1.2) 

<9 t H — curl (u x H) = -curl (v curl H) , divH = 0. (1.3) 

Here x € T d . a torus in R d , d = 2 or 3, t > 0, the unknowns p denotes the density, u = 
(ui, . . . , Ud) G R d the velocity, and H = (Hi, . . . , Hd) £ M. d the magnetic field, respectively. The 
constants fx and A are the shear and bulk viscosity coefficients of the flow, respectively, satisfying 
jl > and 2p, + dX > 0; the constant v > is the magnetic diffusivity acting as a magnetic diffusion 
coefficient of the magnetic field. P(p) is the pressure-density function and here we consider the 
case 

P(p)=ap\ (1.4) 

where a > and 7 > 1 are constants. 

The well-posedness of the Cauchy problem and initial boundary value problems for Ql.ip - Ql.Sp 
has been investigated recently. The global existence of weak solutions to the compressible MHD 
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equations with general initial data was obtained by Hu and Wang [9l|T0] (also see |6| on "varia- 
tional solutions"). From the physical point of view, one can formally derive the incompressible 
models from the compressible ones when the Mach number goes to zero and the density becomes 
almost constant. Based on this observation, Hu and Wang [TT] proved the convergence of the 
weak solutions of the compressible MHD equations (|l.lj) - (|1.3j) to a weak solution of the viscous 
incompressible MHD equations. Jiang, Ju and Li [12] obtained the convergence towards the strong 
solution of the ideal incompressible MHD equations in the whole space by using the dispersion 
property of the wave equation if both the shear viscosity and the magnetic diffusion coefficients go 
to zero. 

In this paper, we shall extend the results on the Cauchy problem in [12] to the periodic case. 
First, we consider the well-prepared initial data for which the oscillations will never appear. We 
will rigorously show the weak solutions of the compressible MHD equations converge to the strong 
solution of the ideal incompressible MHD equations in the periodic domain if both the shear 
viscosity and the magnetic diffusion coefficients go to zero, as well as to the strong solution of the 
viscous incompressible MHD equations. Furthermore, we shall also give the rates of convergence 
which are not obtained in [TT]II2] . Secondly, we consider the case of general initial data. For this 
case the oscillations (acoustic waves) will appear. Comparing with [12] where the Cauchy problem 
was dealt with, the acoustic waves in the current situation will lose the dispersion property and 
will interact each other. Thus, here we have to impose more regular conditions than L 2 on the 
initial data to control the oscillating parts. In addition, we have to assume that the Sobolev norm 
of the oscillating parts is comparable to the magnetic diffusion coefficient in order to deal with 
the general initial data. We will rigorously prove the convergence of the weak solutions of the 
compressible MHD equations to the strong solution of the incompressible MHD equations, as well 
as to the strong solution of the partial viscous incompressible MHD equations. 

To begin our argument, we first give some formal analysis. Formally, by utilizing the identity 

V(|H| 2 ) = 2(H • V)H + 2H x curlH, 
we can rewrite the momentum equation (|1.2|) as 

dt(pu) + div(pu <g> u) + VP = (H • V)H - -V(|H| 2 ) + /EtAu + (fi + A)V(divu). (1.5) 
By the identities 

curl curl H = V divH - AH 

and 

curl (u x H) = u(divH) - H(divu) + (H ■ V)u - (u • V)H, 
together with the constraint divH = 0, the magnetic field equation (|1.3[) can be expressed as 

<9 t H + (divu)H + (u • V)H - (H • V)u = i>AH. (1.6) 

We introduce the scaling 

p(x,t) = p e (x,et), u(x, t) = eu e (x, et), H(x, t) = eR e (x, et) 
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and assume that the viscosity coefficients p, £, and v are small constants and scaled like 

fl = ep e , A = eA £ , i> = ei/ £ , (1.7) 

where e G (0, 1) is a small parameter and the normalized coefficients p e , A £ , and v t satisfy p e > 0, 
2p e + dA £ > 0, and v e > 0. 

With the preceding scalings and the pressure function (|1.4[) . the compressible MHD equations 
([LI) . (fl~5[) . and (fl~6| take the form 

+ div(p £ u £ ) = 0, (1.8) 

d t (p £ u £ ) + div(p £ u £ <g> u £ ) + aV( f )7 = (H £ • V)H £ - iv(|H £ | 2 ) 

+ /i £ Au £ + (// + A £ )V(divu £ ), (1.9) 
3 t H £ + (divu £ )H £ + (u £ ■ V)H £ - (H £ ■ V)u £ = v e AH f , divH £ = 0. (1.10) 

Moreover, by replacing e by y/a/y e, we can always assume a = 1/7. 

Now, we investigate the incompressible limit of the compressible MHD equations (|1.8p - (|1.10[) . 
Formally let e — > in the equations (|1.8p - ([1.10[) . then we obtain from the momentum equation 
(11.91) that p e converges to some function p(t) > 0. If we further assume that the initial datum p 1 ^ 
is of order 1 + O(e) (this can be guaranteed by the initial energy bound (|2.4|) below), then we can 
expect that p = 1. Thus, the continuity equation (|1.8[) gives divu = 0. Furthermore, using the 
assumption 

/i £ ^0, i/->0 as e^0, (1.11) 

we obtain the following ideal incompressible MHD equations (suppose that the limits u £ — > u and 
H £ -)• H exist) 

8 t u + (u • V)u - (H • V)H + Vp + iv(|H| 2 ) = 0, (1.12) 
c*,H + (u-V)H-(H-V)u = 0, (1.13) 
divu = 0, divH = 0. (1.14) 

In Section 3 we shall rigorously prove that the weak solutions of the compressible MHD equations 
(|1.8I) - (|1.10[) converge to, as e — > 0, the strong solution of the ideal incompressible MHD equations 
(|1.12p - (ll,14p for the well-prepared initial data in the time interval where the strong solution of 
()1.12p - (|1.14p exists. Furthermore, the convergence rates are obtained. To show these results, since 
the viscosity coefficients go to zero, we lose the spatial compactness property of the velocity and 
the magnetic field, and the arguments in [11) do not work here. To overcome such difficulty, we 
shall carefully exploit the energy arguments. 

Next, if we assume that the shear and the bulk viscosity coefficients and the magnetic diffusivity 
coefficient satisfy 

/i £ -> p > 0, A £ -> A, v c -> v > as e -> 0, (1-15) 
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then the compressible MHD equations (|1.12|) - (|1.14|l formally converges to the incompressible MHD 
equations (suppose that the limits u e — > u and H e — > H exist) 

d t vL + (u • V)u - fiAu + Vp - (H ■ V)H + iv(|H| 2 ) = 0, (1.16) 

9 t H + (u • V)H - (H • V)u - i/AH = 0, (1.17) 

divu = 0, divH = 0. (1.18) 

In Sections 3 and 4 we shall prove the convergence to the strong solution of the incompress- 
ible viscous MHD equations (|1.16p - (|1.18p for both the well-prepared and the general initial data. 
Furthermore, the convergence rates are also obtained for the well-prepared initial data. For the 
general initial data, we shall also show the convergence to the strong solution of the partial viscous 
incompressible MHD equations (that is, /j, = and v > in (|1.16[) - (|1.18l) ). 

There are a lot of studies on the compressible MHD equations in the literature. Besides the 
aforementioned results, the interested reader can see |131ll~7] on the global smooth solutions with 
small initial data and see 0(211 on the local strong solution with general initial data. We also 
mention the work |25| where a MHD model describing the screw pinch problem in plasma physics 
was discussed and the global existence of weak solutions with symmetry was obtained. 

Before ending the introduction, we give the notation used throughout the current paper. We 
denote the space L|(T d ) by 

L\{T d ) = {/ G L [oc (T d ) : /1 { |/|>i/ 2} G L\ /1 { |/|<i/ 2 } G L 2 }. 

We use the letters C and Ct to denote various positive constants independent of e, but Ct may 
depend on T. For convenience, we denote by H r = H r (T d ) (r G K) the standard Sobolcv space. 
For any vector field v, we denote by Pv and Qv the divergence-free part and the gradient part of 
v, respectively. Namely, Qv = VA _1 (divv) and Pv = v — Qv. 

We state our main results in Section 2 and present the proofs for the well-prepared case in 
Section 3 and the ill-prepared case in Section 4, respectively. 

2. MAIN RESULTS 

We first recall the local existence of strong solutions to the ideal incompressible MHD equations 
(fTTT2]) - (prT4|) in the torus T d . The proof can be found in [U[23]. 

Proposition 2.1 ( [4"ll23|). Assume that the initial data (u,H)| t= o = (uo,Ho) satisfy Uo,Ho G H s 
(s > d/2 + I), and divuo = 0, divHo = 0. Then, there exist a T* G (0,oo) and a unique solution 
(u,H) G L°°([0,r*),P s ) to the ideal incompressible MHD equations (|1.12|) - (I1.14|) satisfying, for 
anyO<T< T* , div u = 0, divH = 0, and 

sup {||(u > H)(t)|| H . + ||(a t u,ftH)(t)|| H .-i + ||VK*)||h.-i} < C T . (2.1) 

0<t<T 

Remark 2.1. The local existence of strong solutions to the incompressible viscous MHD equations 
(fTTT6jl - (pTT8jl was also established in [41123]. 
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We prescribe the initial conditions to the compressible MHD equations (|1.8p - (|1.10j) as 
p e |t=o = Po(x), p c u e | t=0 = p e (x)u e (x) = m e {x), H e | t=0 = H^(x), 
and assume that 

Po > 0, Po € L 7 , Poluol 2 G L\ H G i 2 , divH = 0, m = for a.e. p e = 0. 
Moreover, we assume that the initial data also satisfy the following uniform bound 



l u ol 2 + 2 |H « |2 + ^TTi) ( W - 1 " 7(p 1)) 



dx < C. 



(2.2) 
(2.3) 
(2.4) 



The initial energy inequality (|2.4[) implies that /Og is of order 1 + 0(e). 

Under the above assumptions, it was proved in [5] that the compressible MHD equations (|1.8I) - 
(II.IOP with initial data (|2.2p ~ (|2.4l) has a global weak solution. More precisely, we have 



Proposition 2.2 ( [9]). Let"/ > d/2. Suppose that the initial data (pg,u ,Hg) satisfy the assump- 
tions (|2.3p and (|2.4p . TTien i/ie compressible MHD equations (|1.8p - (|1.10p iwi/i £/ie initial data (|2.2p 
enjoy at least one global weak solution (p e ,u e ,H e ) satisfying 

(1) p e € L°°(0,oo;LT) n C([0,oo),i7) /or aZZ 1 < r < 7, p e |u e | 2 e L°°(0, 00; L 1 ), W G 
L°°(0, oo;L 2 ), and u £ G L 2 (0, T; ff 1 ), ^u 6 G C([0,T],L^ fc ), H e G L 2 (0, T; tf 1 ) n 

C([0,T],i^)/ O r Te (0,oo); 

(2) i/ie energy inequality 

£ e {t) + f V e {s)ds <£ e (0) 
Jo 

holds with the finite total energy 

and the dissipation energy 

V\t) = f [/i e |Vu e | 2 + (M e + A e )|divu e | 2 + i/ e |VH e | a ](t); 

(3) i/ie continuity equation is satisfied in the sense of renormalized solutions, i.e., 

d t b(p e ) + div(6(p e )u e ) + (b'(p e )p e - 6(p e ))divu £ = 

for any b G C (T) smc/i £/ia£ o'(z) is constant for z large enough; 

(4) tte equations $Tty - $FM hold in V (T d x (0,oo)). 

The main results of this paper can be stated as follows. 

Theorem 2.3. Let s > d/2 + 2 and // + A £ > 0. Suppose that the initial data (pg, Ug, Hg) satisfy 
the conditions presented in Proposition Assume further that 

(2.9) 
(2.10) 



(2.5) 
(2.6) 

(2.7) 
(2.8) 



f |po-l| 2 l (K) _ 1 |< 5) dx-+ f |p§ - l\n (H _ 1>g) dx < Ce 2 , 
\\\/ph u o ~ u 1 1 ^2 (T d } < Ce, ||H - H || 2 2(Tti) < Ce 
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for any S £ (0,1), where uq and Ho are defined in Proposition \2.1\ We assume that the shear 
viscosity fj, e and the magnetic diffusion coefficient v e satisfy 

p e = e a , v e = e 13 (2.11) 

for some constants a, ft > satisfying < a + ft < 2. Let (u, H) be the smooth solution to the 
ideal incompressible MHD equations (|1.12[) - (|1.14[) defined on [0, T*) with (u,H)| t= o = (uo,Ho). 
Then, for any < T < T* , the global weak solution (p e , u e , H e ) of the compressible MHD equations 
(|1.8[) - (|1.10p established in Provosition \2. 2\ satisfies 

P e - l| 2 l(|^_n< 5) dx+ f \p*-l\n (lp e_ 1>s) dx<C T e 2 , (2.12) 
||V^u e -u||2 <C T e CT , ||H £ -H||2 T <C T e CT (2.13) 



lL 2 (T rf ) ^ u rt , ||" — j.j.|| L 2( T 
for any t £ [0, T], where o = min{a, ft, 1 — (a + ft)/%} 



The proof of Theorem 12.31 is based on the combination of the modulated energy method, moti- 
vated by Brenier [T] ; the weak convergence method and the refined energy analysis. Masmoudi |19j 
made use of such idea to study the incompressible, inviscid limit of the compressible Navicr- Stokes 
equations in both the whole space and the torus. Comparing with the proof in [19] . here we have 
to overcome the difficulties caused by the strong coupling of the hydrodynamic motion and the 
magnetic field. 

Furthermore, we can use an idea similar to that described above to obtain the convergence 
of the compressible MHD equations (|1.8p - ([1.10[) to the incompressible viscous MHD equations 
(|1.16p - (|1.18p . In fact, we have the following result. 

Theorem 2.4. Let s > d/2 + 2 and /i e + X e > 0. Suppose that the initial data (p^, Uq, Hq) satisfy 
the conditions presented in Proposition \2. 2\ Assume further that 



\pl-lW\ Pl -i\< S) dx+ \pl - ipi (|p ._ 1|>4) dx < Ce 2 , (2.14) 



HVpguS-Uolli^) <Ce, ||H -H ||£ 2(Td) < Ce (2.15) 

for any 5 £ (0, 1) and for some Uq, H £ H s (T d ), satisfying divu = 0,divH o = 0. We also assume 
that the shear viscosity /i e and the magnetic diffusion coefficient v> e satisfy f|l . 15[) . Let (u, H) be 
the smooth solution to the incompressible MHD equations (|1.16l) - (|1.18p with (u,H)| t= o = (uo,Ho). 
Then, for any < T < T** ( T** is the maximal time of existence for (|1.16p - (|1.18p ). the global 
weak solution (p e , u e , H £ ) of the compressible MHD equations (|1.8P - (|1.10p established in Proposition 
\2. 2\ satisfies that Vu e and VH e converge strongly to Vu and VH in L 2 (0,T;L 2 (T d )), respectively. 
Moreover, for any t £ [0,T], we have 

\p t -l\ 2 l [ \ P ^ 1 \<8 ) dx+ ( \p*-l\n (lp *_ 1>s) dx<C T e 2 , (2.16) 



^ - n\\ 2 L2(Td) + \\W- H||| 2(Td) < C T ^= + C T ( l -^—M + ^=^) . (2.17) 

V l*- v v v A* 
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To show Theorem 12 .41 besides the techniques mentioned above, we have to employ a new tech- 
nique, that is, to modulate both the total energy and the partial dissipative energy simultaneously. 
Moreover, the dissipative effect of the viscous terms is also carefully exploited to obtain the desired 
results. 

Remark 2.2. Comparing with Theorem 12.31 we have gotten the better convergence rates (|2.17ll 
than (|2.13|) when the shear viscosity and the magnetic diffusion coefficient tend to some positive 
constants. 

Some results in Theorem l2.4l can be extended to the case of general initial data. More precisely, 
we shall obtain the convergence of the compressible MHD equations (|1 .8[) - (| 1 . 10[) to the incom- 
pressible MHD equations (|1.16l) - (ll.f 81) for the general initial data under the conditions that the 
oscillating parts of the initial data have higher regularity and the Sobolev norm of the oscillation 
parts is comparable to the magnetic diffusion coefficient. This implies that the influence of oscil- 
lations on the magnetic field can be balanced by the diffusive effect of the magnetic field, which is 
one of the new ingredients in our paper. 

To describe the result, we write p e = 1 + eip e and denote 



We will use the above approximation H € (x, t) instead of <p e , since we can not obtain any bound for 
ip e in L°°(0, T; I?) directly if 7 < 2. 

Theorem 2.5. Let s > 2 + d/2 and 2/i + dX > 0. Suppose that the initial data (pq,Uq,Hq) 
satisfy the conditions presented in Proposition \2.2i Moreover, we assume that ^//^Uq converges 
strongly in L 2 to some uo satisfying Quo G H s , Hq converges strongly in L 2 to some Ho with 
J Jd Ho(x)dx = 0, n e | t= o = TLq converges strongly in L 2 to some ipo € H" , and 



for some constant cq > 0. Let (u, H) be the smooth solution to the incompressible MHD equations 
(fTTB|) - (fl~TB|) with (u,H)| t=0 = (u ,H ) G H s (T d ) satisfying u = Pu and divH = 0. Then, for 
any < T < T** ( T** is the maximal time of existence for - (|1.18|) ), the global weak solution 

(/O c ,u e ,H e ) of the compressible MHD equations (|1.8[) - (|1.10[) established in Provosition \2.2\ satisfies 

(1) p e converges strongly to 1 in C([0, T], L\ (T d )); 

(2) VH £ converges strongly to VH in L 2 (0,T; L 2 (T d )); 

(3) H e converges strongly to H in L°°(0, T; L 2 (T d )); 

(4) P(y/fu e ) converges strongly to u in L°°(0,T;L 2 (T d )); 

(5) yfffvf converges weakly to u in iJ~ 1 (0, T; L 2 (T d )). 




By slightly modifying the proof of Theorem 12 .51 we can obtain the convergence of compressible 
MHD equations to the partial viscous incompressible MHD equations when the shear viscosity 
goes to zero and the magnetic diffusion coefficient goes to a positive constant. The partial viscous 




(2.18) 
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incompressible MHD equations correspond to the case of turbulent flow with very high Reynolds 
number (where the viscosity of flow can be ignored, see |f 6j). 

Theorem 2.6. Let s > 2 + d/2. Suppose that the conditions in Theorem \2.5\ hold. Moreover, we 
assume that 



and // = e a for some constant < a < 1. Let (u, H) be the smooth solution to the following 
partially viscous incompressible MHD equations 



d t tt + (u ■ V)H - (H • V)u - ^AH = 0, 
divu = 0, divH = 0, 

with (u,H)| t= o = (uo,Ho) G H S (T 3 ) satisfying Uq = Puq and divHo = 0. Then, for any < T < 
T** (T** is the maximal time of existence for (jf .f 6|) - (|1.1 81) ), the global weak solution (p e ,u e ,H e ) 
of the compressible MHD equations fjl.8[) - fjl.l0[) established in Proposition \2.2\ satisfies 

(1) p e converges strongly to 1 in C([0, T], L\ (T d )); 

(2) VH £ converges strongly to VH in L 2 (0, T; L 2 (T d )); 

(3) H e converges strongly to H in L°°(0, T; L 2 (T d )); 

(4) P(y/fu") converges strongly to u in L°°(0, T; L 2 (T d )); 




(5) y/ffu e converges weakly to u in H 1 (0, T; L 2 (T d )). 

Remark 2.3. The assumption that IIq converges strongly in L 2 to some ipo in fact implies that </?q 
converges strongly to ipo in Lj- 

Remark 2.4. When taking H e = in (|f .1 p ~ (|f .31) . the MHD equations reduce to the classical 
compressible Navier-Stokcs equations. The low Mach number limit problem of the compressible 
Navier-Stokes equations has been investigated extensively, for instance, see [U [31 [51 IS]- The 
interested reader can refer to the survey article [20) for more related results. 

Remark 2.5. We point out that our arguments in the present paper can be applied to the case 
of H e = 0. In this case, we obtain the convergence of the compressible Navier-Stokes equtions to 
the incompressible Euler or Navier-Stokes equations with general initial data, extending thus the 
results in pHT9j . 



In this section, we shall prove our convergence results for the case of well-prepared initial data 
by combining the modulated energy method, the weak convergence method, and the refined energy 



v e -)• v > 0, 2^ + A £ ->• 29 > as e -> 0, 



d t u + (u • V)u + Vp - (H • V)H + -V(|H| 2 ) 




3. Proof of Theorems 12.31 and 12.41 



analysis. 



INCOMPRESSIBLE LIMIT OF MHD EQUATIONS 



9 



Proof of Theorem \2.3l We divide the proof into several steps. 
Step 1: Basic energy estimates and compact arguments. 

By the assumptions on the initial data we obtain, from the energy inequality (|2.5|) . that the 
total energy £ e (t) has a uniform upper bound for a.e. t £ [0,T], T > 0. This uniform bound 
implies that p e \u e \ 2 and ({p e )^ - 1 - j{p e - 1)) /e 2 are bounded in L°°(0, T; L 1 ) and H e is bounded 
in L°°(0,T; L 2 ). Using the analysis in [18], we obtain 

jfd e - jjd e 

which implies (|2 . 12[) and 

p e ^l strongly in C([0, T], L^(T d )). (3.2) 

From the results in [T5], we know that 1 1 u - e 1 1 if — ^ + C|| 'Vu e || ^,2^2 • Furthermore, the fact 
that /9 e |u e | 2 and |H e | 2 are bounded in L°°(0,T; L 1 ) implies the following convergence (up to the 
extraction of a subsequence e n ): 

y/p T u e converges weakly-* to some J inL°°(0, T; L 2 (T d )), 
H e converges weakly-* to someK in L°°(0, T; L 2 (T d )). 

Thus, to finish our proof, we need to show that J = u and K = H in some sense and the 
inequalities (|2.13|) hold, where (u, H) is the strong solution to the ideal incompressible MHD 
equations (fT7T2|) - ([TjH|) . 

Step 2: The modulated energy functional and the uniform estimates. 

We first recall the energy inequality of the compressible MHD equations (|1.8I) - (I1.10I) . i.e., for 
almost all t, there holds 



If \f(t)\ur(t) + m 2 (t) + (Ti e (t)) 2 } +// / / ivuf 

L Jf d L J Jo JT d 

+ (/i e + A e ) f f |divu e | 2 + z/ f f |VH £ | 2 

Jo Jr d jo Jr d 

p |u | 2 + |H | 2 + (n ) 2 l. (3.3) 



1 

< - 

2 Jfd 

The conservation of energy for the ideal incompressible MHD equations (|1.12|) - (|1.14|) reads 



\f [|u| 2 (i) + |H| 2 (t)] =\ f [|uo| 2 + |H | 2 ]. (3.4) 

Z Jfd L Jfd 



2 

Using u to test the momentum equation (|1.9[) , we obtain 



/ (p^-u)(t)+f f p^- [(u.V)u-(H.V)H + Vp+iv(|H| 2 )] 

JT d Jo JT d ^ 

-f f [(p e u e g> u e ) • Vu + (H e • V)H e • u - //Vu 6 • Vu] = / p u • u . (3.5) 

JO JT d JT d 
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Similarly, using H to test the magnetic field equation ([1.10)1 . one gets 

/ (H £ -H)(t)+ f f H e ■ [(u • V)H - (H ■ V)u] + v € f f VH £ VH 

JT d JO JT d JO JT d 

+ 11 [(divu e )H e + (u e • V)H e - (H e ■ V)u e ] ■ H = f H £ H . (3.6) 

Jo Jr d Jr d 

Summing up ([3.3)1 and ()3.4j) . and inserting ()3.5|) and ([3.6)1 into the resulting inequality, we can 
deduce the following inequality by a straightforward computation 



\ J - u| 2 (t) + |H £ - H| 2 (t) + (IT )\t)\ 



+ M £ / / IVu^ + ^ + A 6 ) / / |divu £ | 2 + i// / |VH e | 2 

Jo Jr d Jo Jr d Jo Jf d 

<^f f Vu £ .Vu + / f f VH e .VH - / / p e u £ • [(H ■ V)H)] 
Jo Jt 1 Jo Jj d Jo Ji d 

-ff (H e -V)H e -u+ f f H e • [(u • V)H - (H • V)u] 

JO JT d JO JT d 

+ 11 [(divu £ )H e + (u e -V)H e -(H e -V)u e ] -H+i / f p e u e ■ V(|H| 2 ) 

JO JT d 2 J J T d 

+ 11 P e u e • [(u- V)u + Vp] - / / (p £ u £ 8u £ )-Vu 

JO JT d JO JT d 

+ / [Wf - l)^u e • u] (f) - f [(V^ - l)V^u e • u] (0) 

JT d JT d 

+ \J jd {lV^u e - u| 2 (0) + |H £ - H| 2 (0) + (H^) 2 }. (3.7) 

We first deal with the right-hand side of the inequality ([3.7)1 . Denoting w e = \fffu e — u and 
Z e = H e — H, integrating by parts, and using the fact that divH e = 0, divu = and divH = 0, 
we find that 



/ f ■ [(H • V)H)] - / / (H £ • V)H e • u 

JO JT d JO JT d 

f f H £ .[(u.V)H-(H.V)u]+i/7 p £ u e .V(|H| 2 ) 

Jo JT d L JO JT d 

[(divu £ )H e + (u e • V)H e - (H e • V)u e ] • H 

Iq JT d 

f f p e u £ ■ [(H • V)H)] + f f (H e ■ V)u • H c 

Jo Jr d Jo JT d 

f f (u-V)H-H e - f f (H-V)u-H £ 

Jo Jr d Jo JT d 

(u £ -V)H-H E + f f (H e -V)H-u £ + i / / p e u £ -V(|H| 2 ) 

1 Jo JT d 

/'/ ' 

'0 JT d 



10 JT d Jo JT d 2 Jq J T d 

(1 - p e K • m ■ V)H)] + f f [(H e - H) • V]u • (H e - H) 

'■ JO JT d 
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Jf d 



[(H £ - H) • V]H • (u e - u) - / / [(u e - u) • V]H • (H £ - H) 



JT d 



and 



+ \ f ! (P e -1KV(|H| 2 ) 

< f [ (l-p>M(H-V)H)]+ f W{ S )\\lA\Vu{s)\\ L ~d s 
Jo JT d Jo 

+ Al|w e ( s )||| 2 + ||Z e ( s )||| 2 ]||VH( s )||^d s + f f (Z £ -V)H.[(1-^K] 
Jo Jo JT d 

[ {[(l-V7K]-V}H.Z £ + i f f ( P e -1KV(|H| 2 ) 

IT d z Jo JT d 

f f [p e u e - ((u- V)u + Vp)] - f f (pV®u f )-Vu 

JO JT d JO JT d 

-f f (w f «w f )-Vu+ / f (p e - V^)u e ■ ((u • V)u) 

JO JT d JO Jf d 



[{yfjfu f - - U) • V]U ■ W £ 
/o Jl d 

f [ (^u £ -u).V(M!). 

'0 JT d L 



p e u e ■ Vp 



Jf d 



Substituting (|3.8p and (|3.9p into the inequality (|3.7p . we conclude that 
||w e (0|| 2 L2 + ||Z e (t)||i 2 + ||ff(i)||i 2 +2^ /' / |Vu e | 2 

Jo Jf d 



+ 2(// + A £ ) / / |divu e | 2 + 2j/ 



o JJ d 



|VH e 



o J7 d 



<2C A||w e (r)|| 2 2 + ||Z e (r)|| 2 2 )(||Vu(r)|| L . + ||VH(r)|| L »)dr 
Jo 



||w e (o)|| 2 2 + ||z £ (o)|| 2 2 + ||n^||| 2 



where 



Rl(t)=fi € f f Vu e -Vu + z// f VH £ • VH, 

JO JT d JO JT d 

*§(*)=/ [(V^-i)V^u e -u](*)- / [(^-i)V^u £ -u](o) 

JT d JT d 



JT d 



(p e -^K.((u.V)u), 



i2§(t) = / / (Z e • V)H • [(1 - - /*/ { [(i - • V}H ■ 7f 

Jo Jf d Jo JT d 

R%{t)=f f (l- P e K.[(H.V)H)] + i /V (^-1K.V(|H| 2 ), 

Jo JT d z JO JT d 



(3.8) 



(3.9) 



(3.10) 
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Rl(t)= / / p £ u e -Vp, 

JO JT d 

m)=- f I (V7u e - U ).v(^f ). 



Step 3: Convergence of the modulated energy functional. 

To show the convergence of the modulated energy functional p.lOp and to finish our proof, we 
have to estimate the reminders -Rf(t), i = 1, . . . , 6. 

First, in view of (|3.1[) and the following two elementary inequalities 

\y/x- 1| 2 < M\x- 1| 7 , \x- 1| > 5, 7 > 1, (3.11) 

\y/x- 1| 2 < M\x- 1| 2 , x>0 (3.12) 

for some positive constants M and < 6 < 1, we obtain 



<Me" 



P «-i|<i 



P «-ii>i 



ip--ii<I 



\p e -l\ 2 + M 



Ip 5 -i|>| 



|p e -l| 



(3.13) 



Now, we begin to estimate the terms Rf(t),i = 1, . . . , 6. For the term R\(t), by Young's inequality 
and the regularity of u and H, we have 



\Rl(t)\ < 



i< 



|Vu 



ei2 



VH e | 2 + C T Ai e + C* T ^ 



T d 



(3.14) 



/o JT d 2 J 

For the term i?|(<), by Holder's inequality, the estimate (|3.13[) , the assumption on the initial 
data, the estimate on y/p^u e , and the regularity of u, we infer that 



TO)|<Ce+||u(t)|| L -(/ \Vf-l\ 2 Y([ P>TV 



(u-V)u](t)|| i0 



JT d 



JT d 



<C T e- (3.15) 

For the term i?| (t) , making use of the inequality (|3.13[l , the basic inequality (|2.5[) , the estimates 
on u £ and H e , the regularity of H, the assumption p. lip , and Sobolev's imbedding theorem, we 
get 

\Rl(t)\ <(II[(H • V)H](i)|| L <» + \\vn(t)\\ L ~ ■ \\n(t)\\ L ~) 



10 Jl d 

|VH(t)||i 



JT d 



p*-lH ||u e (r)|Ue||H e (r)|| L 3 



(It 



<C T e(l + (^)-^) + ||Vp ? -l| 



L°°(0,T;L 2 
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<C T e(l + {n t )-*) + CTe{\\ U e ||L2(0 T;i 2) + ||Vu £ || i 2(0 ! T;L2)) 

x (I|H £ ||l2( 0T . L 2) + ||VH e || £ 2 (0:T;i 2 :) ) 
<C T e(l + + C T e[l + • [1 + („«)-*] 

<C T e 1 ~ a/2 + C T e a < C T e a , 



(3.16) 



where a = 1 - (a + £)/2. 

For the term R%(t), one can utilize the inequality (|3.13|) . the estimates on u e and -^/p^u 6 , the 
regularity of H, and p t — 1 = p e — ^fff + yfff — 1 to deduce 

,-t 



|i?|W|<(||[(H.V)H](t)|U» + ||V(|H| 2 )|| £ -)( 



JT d 



'0 JT d 
i -i 



p*\x?\* 



JT d 



<C T e(l + (fx e )~i) < Cre 1 ""/ 2 . 
Using (|2.ip . (|3.ip and (J3T2J) , the term J?§ (t) can be bounded as follows. 

|i*s(*)l= / / P^-Vp 

JO JT d 



{((p'-l)p)(t)-((p'-l)p)(0)}- / (p f -~l)d t p 



JT d 



< 



|P C -1|<^ 



|P £ -1| : 



T 3 



W)\- 



T 3 

T-l 



b(o)l : 



v ^|p«-l|<i 

( 



lp e -ll : 



T d 



\dtp(t)\- 



+ ( ip e -ir)"( / \d tP (t)\^) 

Jo v -V-n>i ' K h d ' 



<\p'--\\> 
<C T (e + e 2/K ) < C T e, 

where K = min{2,7} and we have used the conditions s > 2 + d/2 and 7 > 1. 
Finally, to estimate the term Rl(t), we rewrite it as 



m)=- 



(v^u e -u)-Vi 



u 



JT d 



V7(V^-iK-v(^) 



10 JT d z 



where 



Rk(t) = f I 

JO Jf d 



(3.17) 



(3.18) 



(3.19) 
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fl&w = fj r y - i)ft(^) - / Td [(y - 1) (?))(*) - (y - i)(^))(o) 

Applying arguments similar to those used for and Rl(t), we arrive at the following bound- 

edness 

\Rm\ < \RlM + \Rl2(t)\<C T e. (3.20) 
Inserting the estimates (|3.14|l - (|3.20|l into (|3.10j) and applying Gronwall's inequality, we conclude 

\\^(t)\\l 2 + \\Z^t)\\l 2 + \\U^t)\\h 

<C[||w e (0)||| 2 + ||Z e (0)||| 2 + ||n5||| 2 +C7 T 6-], fora.e. t€[0,T], (3.21) 

where 

C^cxpjc^ [||Vu(r)|| L oo +||VH(r)|| L =o]dr} < +oo. (3.22) 

Now, letting e go to 0, we obtain K = H in L°°(0,T;L 2 ) and J = u in L°°(0,T; L 2 ). The 
inequality (|2.13[) follows from (|2.10j) and (|3.21[) directly. Thus, we complete the proof. □ 



Proof of Theorem \2.4\ For simplicity we assume here that p € = p, A e = A, and u e = v are con- 
stants, independent of e, satisfying p > 0, p + A > 0, and v > 0. The case (jl.7|) can be treated 
similarly. The proof of Theorem 12.41 is similar to that of Theorem 12.31 Since the viscosity is in- 
volved here, we have to modulate the part of dissipation energy in the energy inequality (|2.5|) . We 
state the main different points in the proof here. 

From the basic energy inequality (|2.5j) . we obtain that, for a.e. t G [0,T], p^u^ 2 and ((p e ) 7 — 
1 - l{p e - 1)) /e 2 are bounded in L°°(0, T; L 1 ), H e is bounded in L°°(0, T; L 2 ), Vu e is bounded in 
L 2 (0,T;L 2 ), and VH e is bounded in L 2 (0,T;L 2 ). Therefore, we have 

p f ^l strongly in C([0, T], L 7 2 (T d )), 

and u e is bounded in L 2 (0,T;L 2 ). The boundedness of p e |u e | 2 and |H e | 2 in ^(O^;! 1 ) implies 
the following convergence (up to the extraction of a subsequence e„): 

yfffvf converges weakly-* to some J in L°°(0,T;L 2 {T d )), 
H e converges weakly-* to someK ini°°(0, T; L 2 (T d )). 

Our main task in this section is to show that J = u and K = H in some sense, where (u, H) is the 
strong solution to the viscous incompressible MHD equations (|1.16j) - ()1.18j) . 

Next, we shall also modulate the energy inequality (|2.5p . The conservation of energy for the 
viscous incompressible MHD equations (|1.16[) - (|1.18|) reads 

\( [|u| 2 + |H| 2 ](t) + f j [ M |Vu| 2 + HVH| 2 ] =\ f [|u | 2 + |Ho| 2 ]. (3.23) 
Similarly to Step 2, we use u to test the momentum equation (|1.9p to deduce 



/ (p^ ■ u)(t) + f f P ^ • [(u • V)u - (H • V)H - /iAu + Vp + iv(|H| 2 )] 

JT d JO JT d * 
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/ [(p £ u e 8 u £ ) • Vu + (H f ■ V)H £ ■ u - /iVu £ ■ Vu] = / p e u e ■ u . (3.24) 
Ifd Jfd 

Then, wc test (fTTTu]) by H to infer that 

/ (H £ -H)(t)+ f f H e • [(u- V)H- (H V)u-^AH] + v f f VH £ VH 

JT d JO JT d JO Jj d 

+ 11 [(divu e )H e + (u e • V)H e - (H e • V)u e ] • H = f Hq • H . (3.25) 
Jo Jj d Ji d 

Summing up (O and (1535)1 , and inserting 053)), (JSJEJ with p e = \x and A e = A, ([5341) , and ([535)1 
into the resulting inequality, we deduce the following inequality by a straightforward calculation 

\ {|V7u e - u| 2 (t) + |H £ - H| 2 (t) + (IP ) 2 (t)} 



+ M / / |Vu e -Vu| 2 + ^ 

JO JT d z JO JT d 

ft 



|VH E - VH| 2 + (p + A) / / |divu e | 

JO JT d 



e|2 



+ ^ f [ + ^ I I ivhi 2 

z JO JT d z Jo JT d 

< - T / P e u £ • [(H • V)H)] - f [ (H e • V)H e • u 

Jo JT d JO JT d 



10 JT< 
ft 



10 JT d 



W [(u-V)H-(H-V)H] +- / / p e u £ .V(|H| 2 ) 

^ JO JT d 

[(divu e )H e + (u e • V)H e - (H e • V)u £ ] • H + n f f (1 - p e )u € Au 

Jo Ji d 



+ 11 P e u £ • [(u • V)u + Vp] - / / (p £ u e 8u £ )-Vu 

JO JT d JO JT d 

+ [ [(y/P - l)^U e • U] (t) - / [(V^ - 1)^U £ • U] (0) 

JT d JT d 

+ 5 L { |v7u£ _ u|2(o) + |hc ■ h|2(o) + (n ° )2 l- 

By virtue of (|3.8[) and (|3.9[> . we can rewrite the inequality (|3.26[) as follows 
/7u e -u| 2 (t) + |H e -H| 2 W + (ff) 2 (t)} 



(3.26) 



JT d 



|Vu e - Vu| 5 



+ (M + A) / / |divu £ | 2 + - 

JO JT d z 



2 Jo JT d 
/•t 

v 



|VH £ - VHI 2 



|VH e 



o Jt d 



IVHI 



o Ji d 



-\L { IV7u£ - u|2(o) + |H£ ~ H|2(o) + (n ° )2 ^ 

+ J^(t) + Rl(t) + Rl(t) + R%{t) + R\(t) + Rl(t), 



(3.27) 
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where R\{t) and R\{t), i = 4, 5, 6, are the same as before, and 

RMt) = M f ! (1 - P e )u e A M , 
Jo J7 d 

Rl(t) = f f (Z e • V)H • [(1 - - / / {[(1 - V^)u e ] • V}H • Z e . 

Form the previous arguments on i?|(£) and Rf(t), i = 4, 5, 6, we get 

6 



(3.28) 



Now, we estimate -Ry(i) and Rl(t). Using the inequality (|3.13[) . Holder's inequality, the estimates 
on u c and ^/p^u 6 , the regularity of u, and p e — 1 = p e — yj~ff + ^fff — 1, we obtain 



(3.29) 



\Rm<»\\&<t)\\L~( / / IV7-1| 2 ) 5 (/ / |uf) 5 + (/ / p>f) 

Wo JT d 7 L V J0 JIN 7 V J0 JT d 7 



<C T e(l + -L)<C T ^-. 
V/ 1 



For the term i?|(t), we can make use of (|3.13j) . (|2.5|) . and the estimates on u e and H e , the 
regularity of H, the assumption p. lip , and Sobolev's imbedding theorem to deduce 

<(II[(H • V)H](t)|| L - + ||VH(t)|U» • ||H(t)|| £ c.) 



JT d 

VH(t)|U 
1 



JT d 



•t 

' dr 



||u £ (r)|| i8 ||H e (r)|| L 



II 



<C T e(l + -^) + ||y7^llL-(o,T;^)( J o \\u^r)\\ 2 H1 dT 
<C T e(l + — ) + C T e(||u e || L2( o iT ;i2) + \\Vu e \\ L ^ 0iT . L ^) 



X (l|H e ||i2( 0! T;L 2 ) + l|VH C || i 2( ^T;L 2 )) 

<C T e(l + — ) + C T e[l + /i"'] ■ [l + 
v I 1 

<C T e + C T e/^ <C T e/^- (3.30) 
Now, substituting (|3.28l) - p.30l) into (|3.27[) and applying Gronwall's inequality, we conclude 

w^(t)\\h + \\z e (t)\\h + m e (t)\\h 

<C7[||w £ (0)||| 2 + ||Z e (0)||| 2 + ||n^||| 2 +C7 T e/V7^], fora.e. t e [0,T], 



(3.31) 

where C is defined by (|3~2"2"j) . Combining (|2~T5)) with (|33T1) we obtain (|2~T7|) . Substituting (pPTj) 
into (|X27)1 . we conclude that Vu £ converges to Vu strongly in L 2 (0,T; L 2 (T d )) and VH e to VH 
strongly in L 2 (0, T; L 2 (T d )). This completes the proof of Theorem^ □ 
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4. Proof of Theorem 12.51 



In this section we shall study the incompressible limit of the compressible MHD equations (|1.8I) - 
(11.101) with general initial data. Compared with the case of the well-prepared initial data, the main 
difficulty here is to control the oscillations caused by the initial data. For simplicity, we assume 
here that pf = p, A e = A, and v e = v are constants, independent of e, satisfying p > 0, 2p + d\ > 0, 
and v > 0. 

Proof of Theorem \2.5\ As stated in the proof of Theorem 12.41 we obtain from the basic energy 
inequality ([23]) that, for a.e. t G [0,T], //|u e | 2 and ((p e ) 7 - 1 - j(p e - l))/e 2 are bounded 
in L°°(0, T; L 1 ), H e is bounded in L°°(0,T;L 2 ), Vu £ is bounded in L 2 (0,T;L 2 ), and VH e is 
bounded in L 2 (0,T;L 2 ). Therefore, we have 



and u e is bounded in L 2 (0,T;L 2 ). The fact that p e |u e | 2 and |H e | 2 are bounded in L oo (0,T;L 1 ) 
gives the following convergence (up to the extraction of a subsequence e„): 

\fffu e converges weakly-* to some J in 

H e converges weakly-* to someK ini°°(0, T; L 2 {T d )). 

Our main task in this section is to show that J = u and K = H in some sense, where (u, H) is 
the strong solution to the incompressible viscous MHD equations (|1.16j) - (|1.18j) . The key point is 
to control the oscillations caused by the initial data. This can be done as follows. 

Step 1: Description and cancelation of the oscillations. 

In order to describe the oscillations caused by the initial data, we employ the "filtering" method 
which has been used previously by several authors, see [2l[7l fl8lH9] . 

We project the momentum equation (|1.9|) on the "gradient vector-fields" to find 



p e ^l strongly in C([0, T], L](T d )) 



(4.1) 



d t Q(p £ u e ) + Q[div(p e u e ® u e )] - (2/x + A)Vdivu E + -V(|H £ | 2 ) 




(4.2) 



Noticing p e = 1 + e(p e , we can write (|1.8I) and (I4.2p as 



edtip* + divQ(p £ u e ) = 0, 
ed t Q(p e u e ) + Vip e = eF c 



(4.3) 
(4.4) 



where F e is given by 




+ Q[(H £ .V)H e ]--V(( P T-l-7(p' 



f 



!))■ 



(4.5) 
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Therefore, we introduce the following group defined by £(t) = e rL , r £ R, where L is the operator 
defined on V x {V') d with V = {<f> £ V, J Jd cf>(x)dx = 0}, by 



L 



— div v 
-V<f> 

Then, it is easy to check that e rL is an isometry on each H r x {H r ) d for all r £ R and for all 
r £ R. Denoting 

■r(r) 

we have 



- = - dl vv, ^ = -V^». 



Thus, |4 - A0 = 0. 

In the sequel, we shall denote 



and use the following approximations 



U e = _ . , V £ = £f - 



which satisfy 



Q(v^u £ ) I V eJ V Q(y^u' 



|U £ -U £ || »l 40 as t40. (4.6) 

L°°(0,T;L^+i (!"*)) 



With this notation, we can rewrite the equations (|4. 3[) - (|4.4[) as 



or equivalcntly 



5tU e = -LU C + F* 
e 



<9 t V £ = £( (4.7) 



where (and in what follows) v denotes (0,v) T . 

It is easy to check that F e , given by (|4.5|) . is bounded in L 2 (0,T; i?~ s °(T d )) for some s (so G 
R). Hence, V e is compact in time. Moreover, by virtue of the energy inequality (|2.5p and the 
boundedness of the linear projector P, V e £ L°°(0,T; L-1+ 1 (1 )) uniformly in e. Thus, 

V e converges strongly to some V in L r (0, T; H~ s ' (T d )) (4.8) 

for all s' > s and 1 < r < 00. 

Denote 6 = 2fi + X, £i(t) the first component of £(r), and £2(1") the last d components of 
£(r). If we had sufficient compactness in space, then we could pass the limit in (|4.7[) and obtain 
the following limit system for the oscillating parts 

d t V+Qi(u,V) + Q 2 (V,V)-^AV = 0, (4.9) 
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where u is the strong solution of the viscous incompressible MHD equations (|1.16[) - (|1.18|) . Q\ is a 
linear form of V defined by 

Qi(v,V) = lim - / C(-s)( ° ]ds, (4.10) 

t^ootJq \ div(v®£ 2 (s)V + £ 2 (s)V®v) J 

and Q.2 is a bilinear form of V defined by 

Q 2 (V,V) = lim - [ £(-s)( ° , , )ds (4.11) 

V ' ™r7 v '\ div(£ 2 (s)V®£ 2 (s)V) + ^V(£i(s)V) 2 J V ' 

for any divergence-free vector field v £ L 2 (T d ) d and any V = (0, Vq) T S L 2 (T d ) d+1 . Actually, the 
convergence in (|4.10[) and f|4. can be guaranteed by the following Proposition. 

Proposition 4.1 ( [H]). For all v E L ri (0,T;L 2 ) and V e L r2 (0,T; L 2 ), we have the following 
weak convergences (r\ and r 2 are such that the products are well defined) 

w-\\mC(- t -)( t ° t l=Qi(v,V), (4.12) 

e^o V e y V div(v®£ 2 (f)V + £ 2 (f)V®v) / V ^ 

w-lim£f--H f f ° ]=Q 2 (V,V). (4.13) 

e^o V J V div(£ 2 (i)V®£ 2 (f)V) + ^V(A(f)V) 2 ) 

The viscosity term in the oscillation equations (|4.9[) is obtained by the following Proposition. 



Proposition 4.2 ( [19j). Suppose that the same hypothesis as in Proposition on V holds. 
Then, we have 



?AV= lim - / C(-s)( ° )ds 
2 t->oot 7 ; V 0A£ 2 (s)V / 



(4.14) 



The following propositions, the proof of which can be found in |19j . play an important role in 
our subsequent analysis. 

Proposition 4.3 ( |19|). For all v, V, Vi and V 2 (regular enough to define all the products), we 
have 

f Qi(v,V)V = 0, f Q 2 (V,V)V = 0, (4.15) 
f [Qi(v,V 1 )V 2 + Qi(v,V 2 )Vi]=0 ! (4.16) 

I [e 2 (Vi,V 1 )V 2 +2Qa(Vi,V a )Vi]=0. (4.17) 

Proposition 4.4 ( |19j). Using the symmetry of Q 2 , we can extend the equality (|4.13[) in Propo- 
sition \4^l\ to the case: 

lTo £ ( ~ e) { ( Idiv[£ 2 (|)Vi ® £ 2 (f)V 2 + £ 2 (f )V 2 <g> £ 2 (|)Vi] 
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Moreover, the above identity holds for Vi € £ 9 (0, T; if r ) and V2 € £ p (0, T; H~ r ) with r£l and 
l/p + l/q = 1. iliso, J^. i£| ) can oe extended to the case where we replace V2 in the left-hand side 
by a sequence V| such that V| converges strongly to V2 in L p (0, T; H~ r ). 



Step 2: The modulated energy functional and uniform estimates. 
Let V° be the solution of the following system 



d t V° + fii(u, V°) + Q 2 (V°, V°) - -AV° = 



with initial data 



V u | t=0 = (^,Quo) 



(4.19) 



(4.20) 



where u is the strong solution of the viscous incompressible MHD equations (|1.16p - (|1.18j) with 
initial velocity Uo. From [T5], we know that the Cauchy problem (|4.19[) - (|4.20[) has a unique global 
strong solution. 

In order to prove the convergence results in Theorem 12.51 we have to bound the term 



/~ffu e - u - £ 5 



V 



L 2 (T d ) 



|H £ -H| 



L 2 (T d ) 



2 

L 2 (T d ) ' 



To this end, we first recall the following energy inequality of the compressible MHD equations 



p%t)\u^(t) + \nr(t) + (tf(t)y 



|Vu £ 



JT d 



+ (ji + A) 



f / |divu e | 2 + z/ ( / |VH £ 

Jo JT d JO JT d 



1 

< - 

~ 2 



Poluol 2 



(n ) 2 , for a.c. t G [0,T]. 



(4.21) 



On the other hand, the conservation of energy for the incompressible viscous MHD equations 
<fl~T6|) - (fTTT8|) reads 



If [|u| 2 (t) + |H| 2 (t)]+ / / HVu| 2 + HVH| 2 ] =\ f [|u | 

z JT d JO JT d L JT d 

For the system (|4.19|) . Proposition 14.31 implies that 

/ Qi(u,V )V o = 0, / Q 2 (V°,V°)V o = 0, 

JT d JT d 



|H„I 2 1 



(4.22) 



from which the following equality follows. 

1 
2 



|VV U 2 = 



|V°(^0)| 5 



(4.23) 
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Using £i(-)V° as a test function and noticing p e = 1 + eip e , we obtain the following weak 
formulation of the continuity equation (|1.8[) 

c-l 



A - V> e (t) 



1 



div £ 2 - V° W + div(p £ u e )£i - V° 



VWo- 



-[ r £ 1 ri)ftv°^= 

JO JT d Ve/ 

We use u and £ 2 (|)V° to test the momentum equation (|1.9|) respectively, to deduce 

/ (p e u e -u)(t)+ [ / p £ u e - [(u-V)u-(H-V)H- M Au + Vp+iv(|H| 2 )] 

JT d JO JT d * 

-/ / [(p e u e ® u £ ) • Vu + (H e • V)H e • u - /A7u £ • Vu] = / p e u e • u 

JO JT d JT d 



(4.24) 



(4.25) 



and 



- /' 1 £2 © av ° " (p ' u,) - /' l (p ' u ' 8 " v ( £2 © v ° 

AiVu e ■ v(£ 2 (-V°) + + A)divu e div(£ 2 (-V° 



/O JT d 



JT d 



(H e -V)H e -£ 2 (-)V - 



d -|Hrdiv(£ 2 



Similarly, we test (|l~T0|) by H to get 

ri 



II l + WD v °) - 1 « • ( 4 - 26 ) 



/ (H e -H)(<)+ / / H e -[(u-V)H-(H-V)u-i/AH]+ J / VH f VH 

JT d JO JT d JO JT d 

[(divu e )H e + (u e • V)H e - (H e • V)u e ] • H = f Hq • H . (4.27) 

JT d 



JT d 



Summing up (|4.21l) . (|4.22j) and (|4.23j) , inserting (I4.24p ~ ()4.27[) into the resulting inequality, and 
using the fact 



d t V° • V £ 




we deduce, after a straightforward calculation, the following inequality: 



1 

2 ,/qpd 



V^u £ -u-£ 2 (-W° 2 (t) + |H e -H| 2 (t) 



V u e -u-£ 2 - V 



JT d 



lT-AHV 1 



|VH £ - VHl 



(*)} 



JT d 



±f f (|VHf + |VH| 2 )+( At + A) /7 div(V-u-£ 2 (I)v° 

^ JO JV JO JT d V V e 7 
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~ 2 



/^u £ -u-£ 2 - V 



(0) + |H e -H| 2 (0) 



it-a - v 



! (o)} 



where 
A\{t) 

Am 

Al(t) 
A\{t) 

Am 



Am 



(^-l)V7u e - u + £ 2 - V° (t) 



■fd 



(V^-l)V^u e - u + £ 



V 1 



(0) 



(ir-^)A(-jv u j(t) 



(n £ -^)A(-jv°j (o), 



JT d 



(/ - l)u e Au, 



JP 



2 Jo JT d 



AV° • V e 



Vu e ■ V £ 



JT d 



V 1 



( M + A) 



/ / divu e -div(£ 2 ( 
Jo Jf d ^ v 



|/7 |VV°| 2 + M /7 V£ 2 (1)V° 2 + ( M + A) /V 

z JO JT d JO Jl J ve/ JO JT d 



div£ 2 (-)V° 



/O JT d 

n/ 



p £ u e • [(H • V)H)] - 



JT d 



(H e • V)H e • u 



+ / / H e .[(u-V)H-(H.V)H]+i / / p^.V(|H| 2 ) 

JO JT d L JO JT d 

+ f f [(divu e )H £ + (u e • V)H £ - (H e • V)u e ] • H, 

JO JT d 

■ [ f p £ u e ■ [(u • V)u] - / / (p £ u e ®u e )- (vu + £ 2 (-)v° 

Jo Jf d Jo JT d v v e ' 



7-1 



JT d 



(ff) 2 div(£ 2 - V° , 



JT d 



[Si(u,V°) + Q 2 (V°,V )] -V 



(4.28) 



(4.29) 
(4.30) 

(4.31) 
(4.32) 



(4.33) 

(4.34) 
(4.35) 
(4.36) 



Step 3: Convergence of the modulated energy functional. 

To show the convergence of the modulated energy functional (|4.28|) . we need to estimate the 
remainders Af(t), i = 1, . . . , 8. In the sequel, we will denote by ui e (t) any sequence of time- 
dependent functions which converges to uniformly in t. For convenience, we also denote w £ = 
y/p^u e - u - £ 2 (!)V°, Z £ = H e - H, and * e = IP - L x (\ )V°. 

For the term A\(t), we employ (|2.5[) . (|3.13[) . the regularity of u and (|4.6[) . and follow a procedure 
similar to that for R^if), to obtain 

\A\(t)\ < C T e + uj e (t). (4.37) 
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On the other hand, the term A^it) has the same bound as R%{t) + Rj{t), thus 

|^(i)|<C T e. (4.38) 

To bound the term A\(t), we integrate by parts and use the fact £2(~)V° = V<? e for some 
function q e and Proposition 14.21 to infer 



Jo JT d 



Vu e - V (C 2 - V 



i f [ Au e .£ 2 (-)V° 



cl- T - 



o Jf d 
i 







•V c 



f f AV-V°+ w E (f), 
Jo Jf d 



> + A) 



divu e -div(£ 2 ( -)V 



/O JT^ 



(m + a) 



JT d 



JT d 



AV-V°+w £ (i), 



and 



AV° ■ V e 



2 Jo JT d 



AV-V°+w £ (i). 



2 Jo JT d 

Thus, recalling 6* = 2[i + A, one has 

A|(t)=a; e (t). 
Similarly, it follows from Proposition 14.21 that 

A%{t)=u e {t). 

Recalling (|3.8|) and using Holder's inequality, the inequalities ()2.5[) and (|3.13[) . the regularity of 
H, and Sobolev's imbedding theorem, we conclude 

A\{t)< f [ (l-p e K.[(H.V)H)]+ r||Z e (r)||| 2 ||Vu(r)||^dr 

JO JT d JO 



(4.39) 
(4.40) 



JO JT d 
/•t 



< 



w e (r)||| 2 +||Z e (r)||| 2 ]||VH(r)|| L .dr 
(Z £ -V)H- [(1 - V7)u e + £ 2 (-)v° 

f^J[(i-^K + £ 2 (^)v ] -v}h. 

i f [ (p £ - l)u £ V(|H| 2 ) 

z JO JT d 

' [||w e (r)||| 2 + ||Z e (r)||i 2 ] • [||Vu(r)|| L - + ||VH(r)|| L -]dT 



Z f 



JT d 



(H e • V)H-£ 2 (-)v° - (£ 2 (-)v° • v)h-H 



dr + C T e. 



(4.41) 
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The term A\(t) can be rewritten as 



(w e ®w e ) ■ V(u + £(-)V 



10 JT 
i- 1 



7-1 



|/ ij {^u.«(u + £ 2 (1)v») + („ + £2 (I)v»)»^u'}-v(u + £(I)v») 



10 JT d 
r-l 



u + £ 2 (i V°Wu + £ 2 (! V°)-V(u + £(i V° 



7-1 



/0 JT 

r-l 



i 2 



div(£ 2 f-)v ) - ( 7 - l)£if-)v°n e divf£ 



V 



JT d 



(p< - V7)u £ ■ ((u • V)u) 



(V7u e -u)-V(^ 



(4.42) 



We have to bound all the terms on the right-hand side of (|4.42j) . Keeping in mind that div v = 
and applying Proposition 14. 11 one obtains 



{div(£ 2 (l)vo)^(£ 2 (l)vo) 



7-1 



£l( 7 )v 



2 div(£ 2 (-W° 



7-1, 



o Jf d 
i 



A(l)v°| 2 }.(v° + u) 



c 



o JT* 
i 



div(£ 2 U)V°) ® (£ 2 (?)V°] 



A 7 V 



(V° + u) 



Q 2 (V°, V°) ■ (V° + u) + u e (t) = w e {t). 



>0 JT d 

From Proposition 14.41 we get 



(4.43) 



J J [C 2 Q V° ® V7u e + Vfu e ® £ 2 V° } • V (v + £ 2 (I) V° 

r' 



JT d 



(7-l)£i(- V ffdiv(£ 2 (- V° 



: / 1 ( div K D v ° 8 ^ ue + ^ ue 8 £2 (7) v °) + ^ - ^HD v ° ne ) } 



u + £ 2 f-)V° 



£ 



'0 JT d 

■ (V° + u 



£ 



JT d 



div(£ 2 (f)v°® v^u e + v^u £ ®£ 2 (f)v°) + (7- i)v(£i(f)v°rr 



div(£ 2 (f) V° $ Q(^u £ ) + Q(V^u £ ) <8 ^2(f)v°) + (7 - l)v(£i(^) V°n e 



(v° + s) 



JT d 
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div(£ 2 ) V° <g> P( v / P ? u e ) + P(V^u £ ) <g> £ 2 (7) V 
/ / {2Q 2 (V ,V).(V +Q) + Q 1 (u,V ).(V° + Q)}+ w e (<). 

JO JT d 



(V° + u) 



(4.44) 



'0 JT 

Similarly, we have 
rt 



JT d 



{v^u e ®u + u®^u £ } • V u + £(- V° = / / Q 1 (u,V)V°+o; 6 (t), (4.45) 



JT d 



and 



[ I { - • £ * (7) v ° + £ * (7) v °) 8 u ) ■ v (" + <t) v °) 



JT d 



Q!(u,V u )V u +a; e (t) 



(4.46) 



On the other hand, using the basic energy equality (|2 .5[) and the regularity of u, following the 
same process as for i?g(t), we obtain 



'0 Jt 

The term A\{t) can be rewritten as 



(/-#)u £ .((u.V)u)- / (V^u e -u)-V(^f) <C T e. (4.47) 

JO JT d z 



(4.48) 



[ [ [Qi(u,V°) + Q 2 (V°,V )] -V + cu^t). 

JO JT d 

Substituting ([Q5]) - ([Qg]l into ([3~32"]l . we conclude that 

|i4|(*)| + |A e 7 (t)| < Cy (||w e || 2 + ||^|| 2 )(||Vu|| L oc. + |v£ 2 (^)v | i=c )dT + C T e + u e (t). 

(4.49) 

Thus, we insert the estimates on A\(t) (i = 1, ■ • ■ ,7) into (|4.28[) to obtain 
\f {|w^| 2 + |Zf + \^}(t)+n f f v(u e -u-£ 2 (-W° 

Z J T d J JT d v v e 7 



JT d 



(|VZ £ | 2 + |VH e | 2 + |VH| 2 ) + fa + A) 



JT d 



div(u e -u-£ 2 (-)V° 



< /_ [Hw e (r)|| 2 2 + ||Z e (r)|| 2 2 ] • [||Vu(r)|Uoo + ||VH(t)|U« + ||V£ 2 (-j V°|U~) 

\ j {\^\ 2 + \zr + m 2 }(o)+CT£+^(t)+Ai(t)+Aut), 

Z Jfd 



(It 



where 



Am 



JT d 



(H e • V)H ■ £ 2 1 - ) V — I £ 2 1 - I V • V ) H • H 



(4.50) 



(4.51) 
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Now, to deal with and Ag(t), we denote Hq = J Jd Hq(x)cIx to deduce from the 

magnetic field equation (|1 . 10[) that 



W(x,t)dx 



U c (x)dx = H |T fl 



(4.52) 



The assumption that Hq converges strongly in L 2 to some Hq implies 



U e {x) - n (x)dx 



< 



< T 



d\ 



H e (x) -H (a;) |da; 

U e (x) -H {x)\ 2 dx 



■fd 



as e -> 0, 



whence 



H -> H as e — >• 0. 



(4.53) 



Using Holder's inequality, Sobolev's imbedding theorem, Poincare's inequality, the isometry 
property of C, and (|4.52[) and (|4. 53[) . we find that 



\Ai( t )\+\Am\<- 



IVH 



e 1 1 2 



VU e \\ 2 ](r)dT + u(t) 



k\\<Po\\h + \\Q^\\m) A||VH £ || 2 + ||VH £ || 2 
v Jo 



H1(r)dr. 



(4.54) 



Thus, substituting (|4.54[) into (|4.50|) and using (|2.18|) . we deduce by GronwalPs inequality that, 
for almost all t £ [0,T], 



|w e (t)||l 2 + 



V(t)\\l 2 + \\^(t)\\ 

~, + || Z e^„2 



t- | |n -Mollis +C T e+ sup w £ (s)|, 

0<s<t J 



(4.55) 



where 



C = cxp C 



T . 



|Vu(r) 



|VH(r) 



|V£ : 



<7) 



V°| 



dr > < +oo. 



Step 4- End of the proof of Theorem \2.5\ 

Letting e go to zero in (|4.55[) . we see that H e converges strongly to H in L°°(0, T; L 2 (T d )). 
Hence, K = H. Combining (|4.50[) with (|4.55|) . we can easily prove that VH e converges strongly 
to VH inL 2 (0,T;L 2 (T d )). 

Next, it suffices to prove (4) and (5) in Theorem [231 Noting that P(£ 2 (7)V°) = and the fact 
that the projection operator P is a bounded linear mapping form L 2 to L 2 , we obtain, with the 
help of (|4.55j) . that 

't- 



sup ||P(V^u e )-u|| L 2 = sup P(V7u e -u-£ 2 (-W 

0<t<T 0<t<T \ 



< sup 

0<t<T 



p e u e - u - £ 2 
as e^O. 



V° 



L- 



(4.56) 
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Therefore, (4) is proved. Utilizing (|4.1j) . we deduce from (|f .8)1 that divu e converges weakly to 
in # _1 ((0, T) x T d ). Thus we obtain easily that Qu e converges weakly to in i/ _1 (0, T; L 2 (T d )). 
In view of the fact that u e is bounded in L 2 (0,T; L 2 (T d )) and ^/Jf converges strongly to 1 
in C([0,T],L 2 (T d )), we see that Q(y/fu e ) converges weakly to in i? _1 (0, T; L 2 (T d )). Obvi- 
ously, the fact v^u 6 = P( v / /Fu c ) + Q( y / r ffu e ) implies the weak convergence of \f(fxf to u in 
H- 1 (0,T;L 2 (T d )). The proof of Theorem ED is finished. □ 

Theorem 12.61 can be shown by slightly modifying the proof of Theorem 12.51 and therefore, we 
omit its proof here. 
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